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ABSTRACT 

Context. There is growing evidence that a treatment of binarity amongst OB stars is essential for a full theory of stellar evolution. 
However the binary properties of massive stars - frequency, mass ratio & orbital separation - are still poorly constrained. 
Aims. In order to address this shortcoming we have undertaken a multiepoch spectroscopic study of the stellar population of the young 
massive cluster Westerlund 1. In this paper we present an investigation into the nature of the dusty Wolf-Rayet star and candidate bi- 
nary W239. 

Methods. To accomplish this we have utilised our spectroscopic data in conjunction with multi-year optical and near-lR photometric 
observations in order to search for binary signatures. Comparison of these data to synthetic non-UTE model atmosphere spectra were 
used to derive the fundamental properties of the WC9 primary. 

Results. We found W239 to have an orbital period of only ~5.05 days, maldng it one of the most compact WC binaries yet identified. 
Analysis of the long term near-lR lightcurve reveals a significant flare between 2004-6. We inteipret this as evidence for a third mas- 
sive stellar component in the system in a long period (>6 yr), eccentric orbit, with dust production occuring at periastron leading to the 
flare. The presence of a near-IR excess characteristic of hot (~ 1300 K) dust at every epoch is consistent with the expectation that the 
subset of persistent dust forming WC stars are short (<1 yr) period binaries, although confirmation will require further observations. 
Non-LTE model atmosphere analysis of the spectrum reveals the physical properties of the WC9 component to be fully consistent 
with other Galactic examples. 

Conclusions. The simultaneous presence of both short period Wolf-Rayet binaries and cool hypergiants within Wd 1 provides com- 
pelling evidence for a bifurcation in the post-Main Sequence evolution of massive stars due to binarity. Short period 0+OB binaries 
will evolve directly to the Wolf-Rayet phase, either due to an episode of binary mediated mass loss - likely via case A mass transfer 
or a contact configuration - or via chemically homogenous evolution. Conversely, long period binaries and single stars will instead 
undergo a red loop across the HR diagram via a cool hypergiant phase. Future analysis of the full spectroscopic dataset for Wd 1 will 
constrain the proportion of massive stars experiencing each pathway; hence quantifying the importance of binarity in massive stellar 
evolution up to and beyond supernova and the resultant production of relativistic remnants and X-ray binaries. 

Key words. stars:evolution - stars:early type - stars:binary:general 
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spectroscopic observations of massive stellar aggregates (e.g. 
Cyg OB2, NGC6231 & 30 Dor - Kiminki et al. [20^ Sana et 
al. (2008 i Bosch et al. 2009 ) whereby homogeneous stellar pop- 
ulations may be efficiently sampled. These studies reveal a uni- 
formly high binary fraction for OB stars (>40%), which, when 
corrected for the observational bias against long period systems 
results in a fraction consistent with unity for at least one region 
(Cyg OB2; Kobulnicky & Fryer 2007 ). 

Another compelling target is the young (4-5 Myr) mas- 
sive (~10^ Mq) galactic cluster Westerlund 1 (Wd 1; Clark 
et al. I2005I I. which contains a significant (>>100) popula- 
tion of massive (M,„„„,/ >30 M^; Ritchie et al. I2010al l stars. 
Multiwavelength observations of the cluster revealed a num- 
ber of indirect binary markers amongst the OB supergiants and 
Wolf-Rayets (WR), with a binary fraction for the latter in- 
ferred to be at least 70% (and potentially consistent with unity; 
Crowther et al. f2006al Clark et al. 120081 Dougherty et al. Mil. 
Motivated by these findings we have undertaken a systematic 
radial velocity (RV) spectroscopic survey to identify and char- 
acterise the binary population of Wd 1 ; a detailed description of 
the project methodology, reduction technique and preliminary 
results may be found in Ritchie et al. (2009 ). 

In this paper we present an analysis of the properties of the 
WC9d star W239 (= WR F, WR77n), for which binarity was sug- 
gested by the presence of a hard X-ray spectrum and a pro- 
nounced near-mid IR excess due to warm circumstellar dust 
(Skinner et al. [20061 Crowther et al. 2006ai. Including W239, 7 
of the 8 WCL {- WC8-9) stars within Wdl demonstrate one or 
both of these observational properties. Such phenomena are typ- 
ically attributed to the wind collision zone in a massive binary 
system where material is first shock heated to high temperatures, 
with the concommitant high densities subsequently permitting 
downstream dust condensation upon cooling (e.g. Piilutskii & 
Usov [T976l Alle n et al. [19721 Usov et al. lTWTl Williams & van 
der HuchtTW"). 

We present a brief description of the spectroscopic and pho- 
tometric datasets employed in this work in Sect. 2, as well as a 
full RV solution - to the best of our knowledge only the fifth for a 
dusty WC binary after WR 1 13, 137 & 140 (Massey & Niemela 
[19811 WilHams et al. 1990; Lefevre et al. 2005 ) and th e LMC 
system HD 36402 (=Br31; Moff^at et al. 1990. WilUams [20171 ). 
The results of a tailored non-LTE atmospheric analysis of the 
WR component of the binary is also included in Sect. 2, while 
we discuss the implications of our findings in Sect. 3, before 
summarising them in Sect. 4. 

2. Data Reduction & Analysis 

2.1. Spectroscopy 

Eleven epochs of spectra were obtained between 2008-9 with 
the Fibre Large Array Multi Element Spectrograph (FLAMES; 
Pasquini et al. 2002 ) mounted on VLT UT2, with the GIRAFFE 
spectrograph operated in MEDUSA mode (HR21 setup), result- 
ing in a resolving power of ~ 16200 between 8484-9001 A. Full 
details of data reduction and analysis may be found in Ritchie et 
al. ([2009l l. 

While Clark et al. (2010) found that W239 was not subject 
to secular variability between 2001-9, Ritchie et al. (2009 ) report 
the presence of RV variations in the C m 8500 & 8664 A lines in 
the 2008 spectra (e.g. Fig. 1). As a result RV measurements of 
all 1 1 epochs of observations were cariied out using the Nelder- 
Mead simplex method to fit Lorentzian profiles to the strong C iii 
lines, providing a robust RV determination that is unaff'ected by 



small-scale line-profile variations due to wind effects. The re- 
sults of this process are summarised in Table 1, with sample fits 
presented in Fig. 1. 

An obvious explanation for the pronounced RV shifts is that 
they are due to reflex motion induced by an unseen companion 
to the WC9d primary. A period search of the RV data strongly 
favours a candidate period of 5.05days. An alias of this period 
at 3.04 days may not be excluded from the current data but re- 
sults in larger (O-C) residuals and appears further disfavoured 
because wind-driven mass loss in the WC phase acts to widen 
the orbit, implying that W239 would have had to have passed 
through a particularly compact configuration prior to its current 
state; we return to this point in Sect. 3.3.1. Consequently, we 
adopt a period of 5.05 days for the remainder of this paper and 
emphasise that the analysis explicitly excludes longer orbital pe- 
riods (> 6 days) and that any residual uncertainty in the period 
does not materially affect our conclusions. 

We show the resultant cur-ve derived from a Levenberg- 
Marquardt fit to the RV data in Fig. 2, with an (orbital) period, 
Por/,=5.053+0.002 dayfl and semi-ampHtude Kn'c=39.7+4.8 
kms"' (Table 2). Eccentricity was fixed at zero, due to the short 
period and expectation that binary interaction would have cir- 
cularized the orbit. We note that deviations from the RV fit are 
greater in the 2009 data (~10-15kms"'), suggestive of either a 
small systematic error in the period determination or wind vari- 
ability superimposed on the orbital motion. The systemic veloc- 
ity derived from the Cm 8664A line is -60.5+3.5 kms with 
the C III 8500A line offset redwards as noted by Ritchie et al. 
(2009 ). While slightly higher than the mean derived from 10 su- 
pergiants by Mengel & Tacconi-Garman (2009 ) it is within the 
range found for other cluster members from our full FLAMES 
dataset (Ritchie et aL [20091 [20T0bl l. 

It might, however, be supposed that such line profile vari- 
ability was the result of rotating global wind structures. Indeed, 
periodic spectral variability attributed to such a cause has previ- 
ously been reported in the optical lines of the apparently single 
WRs WRl & 6 (Morel et al. [19971 Chene & St-Louis 2010}. 
However, we consider this unlikely for several reasons. Firstly, 
the RV line profile variability in W239 encompasses the whole 
profile. This contrasts to the behaviour of WR 1 & 6, where 
superposition of additional transient emission components on a 
stationary underlying profile leads to shifts in the line centre of 
gravity or 'skewness'. This is also the case for the apparently 
single WC8 star WR K within Wdl, for which we find line pro- 
file variabiUty in the line core but no global RV shifts (Ritchie 
et al., in prep.). Secondly, the 5.05 day period present in W239 
is coherent over ~83 cycles; in comparison, the quasi-period of 
WR 1, attributed to wind asymmetries, is only stable for <4 cy- 
cles (Chene & St. -Louis 2010). Finally, this wind variability is 
accompanied by photometric modulation with an identical pe- 
riod in both WR 1 & 6, but Bonanos (120071 ) demonstrate that 
such behaviour is absent in W239. 

Moreover, phase resolved near-IR spectral analysis of 5 WR 
binaries reveals that variability due to wind-wind interaction 
driven asymmetries - rather than binary reflex motion - is in any 
event phased on their known orbital periods (Stevens & Howarth 
1999). We are thus confident that the 5.05 day period identified 
in the RV curve of W239 reflects an underlying binary orbital 
period of the system. 

Motivated by this finding we re-examined the FLAMES 
spectra of W239 to search for the signature of a companion(s), 
noting that the shorter wavelength data were of too low S/N and 

' Errors are internal to the fit. 
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Table 1. Journal of observations. 



MJD° 


Phase (Orbits) 


RV(C III 8664) 


0-C 






(kms"')'' 




54646.1846 


0.000 


-47.4 


-1.0 


54665.0356 


0.731 (3.73) 


-31.1 


-0.5 


546/i.lj4j 


0.9jO (4.y4) 


-ly.l 


5.8 


54692.0423 


0.076 (9.08) 


-62.1 


1.4 


54713.0107 


0.226 (13.23) 


-87.9 


4.6 


54724.0818 


0.417 (15.42) 


-89.8 


2.0 


54728.0554 


0.203 (16.20) 


-91.0 


-1.8 


54734.0613 


0.392 (17.39) 


-85.5 


9.3 


54965.1768 


0.132 (63.13) 


-92.6 


-16.6 


54969.3198 


0.952 (63.95) 


-22.2 


15.0 


55063.0575 


0.503 (82.50) 


-87.1 


-10.5 



"Modified Julian day at the midpoint of two 600s integrations 
(2x500s, 04/09/2008; Ix600s+lx700s, 19/09/2008), ''Phase taking 
ro=54646.1846, elapsed orbits for an orbital period of 5.053d, "RV 
listed for the C m 8664A line, with fitting errors < 4km s"' at all epochs. 



Table 2. Orbital parameters 



Parameter 


Value 


P (days) 


5.053 ± 0.002 


Kwc (km s"') 


39.7 + 4.8 


7 (km s"') 


-60.5 + 3.5 




< 30° 


e 


(fixed) 


a sin i (R©) 


3.97 ± 0.48 



resolution to permit such an analysis and the longer wavelength 
data are dominated by a near-IR excess due to hot dust (Sect. 
2.2 & 2.3). The FLAMES spectra cover Pall-16, which are ex- 
pected to be in absorption in late O/early B stars. Unfortunately 
Pall falls on the red wing of a Civ emission line from the 
WC9 primary, while Pal3 and Pal6 fall at the same wavelengths 
as strong Cm lines. Given the decrease in the strength of the 
Paschen series to higher transitions, Pal2 (8750A) appears to 
be the best remaining binary diagnostic. Unfortunately, while 
this transition is clear of any WR emission lines it does overlap 
the interstellar C2 Phillips (2-0) R/Q bands, which are visible 
in our spectra of W239 and other cluster members, complicat- 
ing analysis (cf. Ritchie et al. l2010al ). This region of the spectra 
is presented in Fig. 3, and we note that W239 appears to share 
the weak, broad Pal2 profile found for this transition (and the 
wider Paschen seiies) in the O-nO spectroscopic binary W30a 
(Clark et al. 2008; Negueruela et al. 2010b 1; we therefore re- 
gard this as tentative evidence for a hot massive companion(s). 
Unfortunately, the weakness and breadth of the line, combined 
with the low S/N and presence of interstellar features precludes 
a search for RV variability, although pre-empting Sect. 2.3 & 3. 1 
one might expect any reflex motion in a putative early type com- 
panion to be of low amplitude, such that it would be difficult to 
identify at the resolution of these data. 



Table 3. Summary of near-IR photometry 



Date 


H 


nbK 


19/05/99 


7.61 


_ 


28/06/04 


8.21 


_ 


29/06/04 


_ 


6.78 


19/06/04 


8.18 


- 


24/09/04 


7.97 


7.28 


03/06/05 


7.69 


5.78 


27/07/05 


7.67 


_ 


28/07/05 


7.64 


_ 


29/07/05 


7.63 


_ 


31/07/05 


7.64 


5.08 


12/08/05 


7.73 


_ 


16/07/06 


7.56 


- 


04/08/06 


7.59 


_ 


05/08/06 


7.52 


- 


06/08/06 


7.59 


_ 


08/08/06 


7.53 


- 


09/08/06 


7.60 


- 


10/08/06 


7.55 


- 


31/05/07 


- 


6.34 


05/09/09 


7.67 




06/09/09 


7.67 




18/06/09 


7.63 




30/04/10 


7.77 


6.29 


19/07/10 


7.68 




20/07/10 


7.65 




21/07/10 


7.68 




22/07/10 


7.68 





The errors on H band data of ±0.05mag and on nbK band of ±0. 17mag, 
with the exception of the data from 24/09/04 where they are ±0.05 in 
both bands (Crowther et al. l2006al . 



2.2. Photometry 

Wd 1 has been the subject of a long term broad band I, H & 
K and narrow band K0 continuum observational campaign since 
2004, utilising the 0.6-m and 1 .6-m telescopes at Pico dos Dias 
Observatory (LNA/Brazil). These were equiped with an Hawaii 
1024x1024 pixel camera, which gave a field-of-view of 8x8 ar- 
cmin and 4x4 arcmin respectively. Standard image processing 
techniques were adopted (bias subtraction, flatfielding and sky 
subtraction) piior to data extraction via aperture photometry un- 
der IRAF. Finally, absolute I band magnitudes were obtained via 
calibration to the data presented by Piatti et al. ( 1998 1, while the 
near-IR photometry was calibrated via comparison to 2MASS 
data. 

Unfortunately, W239 was saturated in the broad K band, 
leaving just 5 epochs of observations in the narrowband l.lAjjm 
filter. However 25 individual H band observations between 
2004-10 were viable; both datasets are plotted in Fig. 4 and sum- 
marised in Table 3, along with photometry from Crowther et al. 
(2006a). These are further supplemented with the I band data 
from Piatti et al. (1998) and Clark et al. (2005 1 and 10 epochs 
of new observations. Unfortunately, archival 2MASS data were 
associated with large uncertainties and so were not employed. 

These data provide no evidence for secular variability in the 
I band between 1995-2010. However, dramatic near-IR variabil- 
ity is present, with a significant brightening observed between 
2004-6 (AH~0.6 mag, AnbK~2.2 mag). This was also accompa- 
nied by the system becoming redder, with (H-nbK)~0.69+0.07 
on 2004 May 1 and 2.6+0.2 on 2005 July 31. Subsequently, 



- henceforth nbK, centred at 2.14//m, FWHM=0.023/im. 



4 



J. S. Clark et al.: The WC9d binary W239 



W239 appears to have faded in the nbK band, albeit to a higher 
quiescent flux than found during the initial phase of the outburst 
- with the inevitable corollaries that the data are both sparse and 
are associated with comparatively large uncertainties. In contrast 
the H band luminosity was observed to increase by a smaller 
amount, but appears to have remainded at the peak luminosity 
for the following 4 year period. This was also accompanied by a 
modest reduction to (H-nbK)~ 1.3+0.2 by 2010. Unfortunately, 
lack of contemporaneous data did not permit us to determine if 
the outbursts peaked at the same time in both bands. 

Similar photometric behaviour - a transient outburst charac- 
terised by an increase in IR continuum excess with wavelength - 
is observed in other WC binary systems where dust production 
rates are modulated on the orbital period and peak during peri- 
aston passage (e.g. Williams et al. 119941 120011 van der Hucht 
et al. 120011 ). Consequently, we infer the presence of an addi- 
tional long period tertiary companion to W239, with the lack 
of a second outburst suggesting a period of at least ~6 yr for 
this component. However, the broadband near-IR colours at all 
times are clearly indicative of the continued presence of hot dust 
(e.g. Crowther et al. 2006a); in this respect W239 bears a close 
resemblance to WR48a, which also shows both persistent and 
episodic, enhanced dust emission. WilUams et al. ( 120031) also 
interpret this as evidence for a triple system, with a short pe- 
riod binary (Pmb <1 yr) component responsible for continous 
dust formation and an interaction with a tertiary companion in 
a wide, eccentric orbit (Poit >30 yr; R Williams, priv. comm. 
2010) leading to the outburst in 1980. We return to this in Sect. 
3. 

However, in this regard we highlight that despite the sparse 
sampling, the decay in flux from the outburst peak in the nBK 
band lightcurve appeared to occur more rapidly than that of the 
the H band. In other periodic/episodic dust producing systems 
the opposite is observed, with more extended decays at progres- 
sively longer wavelenths due to the cooling of the dust as it is 
carried from the system by the stellar winds. Unfortuntely, given 
the limitations of our current dataset - in particular the lack of 
longer wavelength data - we cannot presently explain this appar- 
ent difference in behaviour 

Finally, optical observations by Bonanos (2007) between 
2006 June 15-July 25 revealed no evidence for eclipses or elip- 
soidal modulation in the lightcurve. We repeated the experiment 
with a higher cadence - ~ 1500 individual I band observations be- 
tween 2008 July 8-17 - but again found no evidence for modula- 
tion on the orbital period although, as with Bonanos, intra-night 
variability of AI< 0.08mag. over short (~hr) timescales was ob- 
served. 

2.3. The physical properties of W239 

In order to investigate the nature of the WC9 primary of 
W239 we have employed the non-LTE model atmosphere code 
CMFGEN (HilHer & Miller 119981 11999] ). To accomplish this, 
in addition to the data described in Sect. 2.1 we also employed 
the spectra presented by Clark et al. (12002) and Crowther et al. 
(I2006al ) and we refer the reader to those works for a description 
of the observational and reduction procedures employed. While 
these data span a significant period of time, no secular spec- 
tral evolution is apparent (Sect. 2.1) and there is an excellent 
correspondence between regions of the spectrum that have been 
multiply sampled. The methodology adopted follows that of the 
analysis of the WC9 star WR103 by Crowther et al. (2006b). 
Given the significant reddening to W239 some of the UV/optical 
transitions utilised by Crowther et al. were unavailable to us; 




-500 -400 -300 -200 -100 100 200 300 400 500 

velocity (km/s) 

Fig. 1. Top panel: sample C in 8664A emission line profiles for 
W239 clearly indicating the binary induced radial velocity shift 
relative to the systemic velocity (dashed line). Lorentzian profile 
fit to emission line as utilised for the radial velocity determina- 
tion overplotted. Note the presence of an additional C iii emis- 
sion feature in the blue wing at 8653.2 A; this was found to have 
no effect on the RV determination. 




-100 I ' ' ' ' ' ' 1 

-0.2 0.2 0.4 0.6 0.8 1 1.2 

phase {T„ = 54646.18456) 



Fig. 2. Radial velocity curve for W239, phased on the orbital 
period. Velocity determined from the C iii 8664 A line. 

thus our primary diagnostics were C uAA 6678, 7230, 9900 A, 
Cm AA 6740, 8500, 9700, 21100 A and Civ^ 14300 A which, 
when combined with local continuum levels enabled an estima- 
tion of temperature, wind properties and stellar luminosity. 

One complication was the significant dilution of the spec- 
trum at long (> 1 yum) wavelengths by emission from hot dust, 
as well as a potential further contribution from a massive com- 
panion(s). The former was accounted for by the inclusion of a 
simple greybody model for the dust emission, assuming a sin- 
gle temperature of 1300 K. To represent the latter we utilised a 
Kumcz LTE model with T~35 kK - corresponding to an 07 V 
star - and logg ~4. We adopted an equal light ratio at 8500 A 
for the WR and OB companion(s). The relative contribution of 
the 07 V star was determined from the analysis of Martins & 
Plez (2006) who found My —4.7 a nd we adopted (V-I)o —0.47 
(appropriate for an O star; Cox et al. 1200 H which yielded M/ ~- 
4.23, while our WC9 model has My —4.3 and (V-I),, —0.05 
leading to M, —4.25. 
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WR F 




0.8 - Pa12 He I Pall- 



8740 8760 8780 8800 8820 8840 8860 

wavelength (A) 

Fig. 3. Comparison of the 8700-8900 A spectra of W239 (solid 
line) and the O+O SB2 binary W30a (dotted Hne). W239 appears 
to demonstrate a similar weak, broad Pa 12 profile at -8500 A as 
W30a. The weakness of this feature is apparent upon considera- 
tion of the single 09.5 lb cluster member W74. 



1 




1 1 
■ 


E(T3-V) = 3 mag 






D=5kpc 










r 




3-6 


3.8 4.0 


4.2 4.4 



log X [A] 



Fig. 5. Comparison between the broadband photometry of W239 
and the reddened WC-i-O spectral energy distribution exclud- 
ing (green dashed line) and including (red solid line) the simple 
1300 K graybody dust model. Optical photometry from Clark et 
al. (2005; 2001 August 21) and the near-IR from 2005 June 3 & 
July 31 (this work). 



1998 



2006 3008 



I band 
H band 
nbK band 



-a 



6 - 



7 - 



FLAMES 
optical 

III I 



SxlO" 



5.2x1 0* 



5.4x10" 



Modified Julian Date [-2400000. 5j 



Fig. 4. Optical-near IR lightcurve for W239. The locations of 
optical and near-IR spectroscopic, radio and X-ray observations 
discussed in the text are also indicated. 



The results of this model are presented in Table 4. 
Encouragingly, they bear close resemblance to the parameters 
of other WC9 stars analysed in an identical manner (Crowther et 
al. l2002ll2006bl l. The resulting best fit synthetic spectiTim (Figs. 
5 & 6) provides an acceptable match to the spectrum from the 
optical to near-IR (~0.6-2.2 yum) and photometry for V through 
K for a distance, d = 5 kpc and reddening, E(B-V) = 3. The resul- 
tant composite synthetic spectrum demonstrates no pronounced 
spectral features from the putative companion; consistent with 
the current observational data. Unfortunately, W239 is saturated 
and/or blended in both Midcourse Source Experiment (MSX) 
and Spitzer Galactic plane surveys; hence more sophisticated 
modeling of the IR excess, although of considerable interest, 
appears unwarranted at this time. At longer wavelengths, the 
3.6 cm radio flux predicted from the model is a factor ~5 smaller 
than observed (~0.3 mJy; Dougherty et al. 2010), suggesting 
the possibility of an additional contribution resulting from the 
wind collision. Finally, we note that our inability to unambigu- 



Table 4. Summary of physical properties of the WC9 primary of 
W239 



Property 



Value 



T, (kK) 
log(L/Lo) 

log(M/V7)(Moyr-') 
/ 

v„ (kms ') 
Xc 



52 
4.95 

3.8 
-4.5 

0.1 
950 
0.25 



ously assign contributions to the continuum from the primary 
and companion(s) precludes us determining an accurate current 
mass for the WC9 star via the H-free mass-luminosity calibra- 
tion of Schaerer & Maeder (11992) , although we consider 10 Mq 
to be a conservative upper limit. 

3. Discussion 

3.1. The nature of the W239 system 

The combination of cyclic RV variability, an IR excess, hard X- 
ray emission when single WC stars are not known to emit at 
such energies (e.g. Oskinova et al. 1200311 and the tentative spec- 
troscopic detection of the luminous companion(s) provides com- 
pelling evidence for the multiplicity of W239. Indeed, the mis- 
match between the 5.05 day RV periodicity and long term pho- 
tometric behaviour is suggestive of a triple system. But what can 
be deduced about the putative companions from a synthesis of 
these observational datasets? 

Unfortunately, the nature of the binary companion in the 
5.05 day orbit is somewhat uncertain (in the absence of an RV 
curve from the putative Pa 12 absorption line we may not as- 
sign it to either a binary or tertiary component with certainty). 
Nevertheless, assuming a mass for the WC9 component of ~8- 
10 Mq (Sect. 2.3), the combination of the observed RV semi- 
amplitude, Kwc^ lack of eclipses and stellar radius derived from 
modeling imply a lower limit to the mass for the unseen com- 
panion of ~2 Mq (for an inclination of ~45°). However based 
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Fig. 6. Comparison of the rectified synthetic spectrum of W239 
(red line) to observations (black line), acounting for the WC9 
line dilution from the unseen companion and warm dust. Spectra 
were obtained in 2001 June (1.5m/B&C), 2005 June (NTT/Sofl) 
and 2008 July (VLT/FLAMES). Note the two emission features 
shortwards of 6500A are the result of incomplete sky subtrac- 
tion; telluric contamination is also manifest in the observed spec- 
trum beween -7200-7400 A, L35-L50//m and L77-L8jL/m. 

on evolutionary considerations - specifically the requirements 
placed on the system mass ratio for merger to be avoided during 
interaction (Sect. 3.3) - we suspect the actual companion mass 
is significantly larger, with the consequent constraint that the bi- 
nary is observed at low inclination. 

As an illustrative example, an hypothetical 10 (WC9) + 
20 - 30 Mq (BO - 07 V) binary observed at -30 - 10° would 
replicate the observed kms"'. A putative 07 V com- 

panion would be consistent with our assumption that the WC9 
primary and companion contribute equally at 8500 A, while a 
fainter BO V star would also satisfy this condition while also ac- 
comodating the presence of a tertiary companion of comparable 
spectral type. The resultant orbital separation of only <0.2 AU 
would then be an order of magnitude smaller than that of the 
archetypal dust producer WR140 (Porb ~7.94 yr) during perias- 
tron. We emphasise that neither scenario represents an attempt 
to provide a fully self consistent model for the W239 system; a 
goal which must await future, stricter observational constraints. 

Would such a putative binary be consistent with the cur- 
rent X-ray data for W239? In Table 5 we summarise the prop- 
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erties of known WC binary systems displaying phenomena as- 
sociated with wind collision zones (X-ray and/or dust produc- 
tion). Unfortunately, no counterpart to an X-ray bright WC9 bi- 
nary has yet been classified, although the WC8-I-07.5 III binary 
WRl 10 suggests that companions of comparatively late spectral 
type may yield detectable X-ray emission in conjunction with 
WCL stars, although we note that this system in not a dust pro- 
ducer. 

With Lx ~ 2 X 10^^ ergs^', W239 is an order of magnitude 
fainter than WRl 1 (Clark et al. 120081 1. but is also a much more 
compact system. For the binary configuration described above, 
the orbital separation is only an order of magnitude larger than 
the stellar- radii of the WC9 and 07 V stai's (Martins et al. 20051 1; 
hence one would expect that neither wind will have reached ter- 
minal velocity at the position of the shock, likely reducing the 
resultant X-ray luminosity in comparison to WRll. Moreover, 
one might also expect radiative braking of the WR wind by the 
O star (cf. V444 Cyg; Owocki & Gayley[T995]l, further suppress- 
ing X-ray emission. 

The X-ray observations of W239 were made in two blocks, 
one of 18 ks (-4% of orbit) on 2005 May 22 and 42 ks (-10% 
of orbit) on 2005 June 20, a separation of 5.8 cycles (Clark et 
al. 12)08). No variability in flux was found within or between ei- 
ther observation. While X-ray variability as a function of orbital 
phase has been observed in WRll (e.g. Schild et al. 2004) that 
system is observed at higher inclination, such that significant or- 
bital variation in the wind properties and hence line of sight ex- 
tinction are present, which would be absent for the putative low 
inclination configuration proposed here. Consequently, one can 
at least qualitatively understand the X-ray properties of W239, 
although quantitative confirmation would require a tailored hy- 
drodynamical simulation, which currently appears unwarranted 
given the lack of observational constraints. 

Finally, we briefly turn to the nature of the putative 
Porb >6 yr tertiary component suggested by the IR lightcurve. 
Unfortunately, the lack of contemporaneous spectroscopy over 
this period precludes any firm conclusions, save that a <07 V or 
~09 III star are suggested by the twin constraints imposed by the 
age of Wdl and the requirement of the star to support a substan- 
tive wind to permit dust production, noting that we might also 
expect such a star to contribute to the observed X-ray emission. 

3.2. Dust production in close binaries 

Our current IR observations are consistent with the hypothesis 
that hot dust is always present within the W239 system. Monnier 
et al. (J999) suggested that the subset of the persistently dusty 
WC stars were short orbital period ( < 1 yr) binaries, where dust 
formation occurs in long lived Archimedean spirals that result 
from high density material produced at the wind collision zone 
being carried outwards by a uniformly expanding wind until 
temperatures are low enough to permit condensation (Tuthill et 
al., 1999). With a ~5.05 day orbital period - comparable to that 
of HD 36402 in the LMC (Moff^at et al. 1990, Williams ISOTlT) 
and a factor of ~6 smaller than WRllifl the next most com- 



3 The physical properties of the WC8 stai- in WRll - T, =57.1 kK, 
log(L/Lo)=5.0, log( dotMl ^)=-A.5'i, v„„ = 1550 kms"' (De Marco et 
al. l2000t - are similar to those of W239, save for a faster wind. In con- 
trast, earlier WC stars are significant hotter, more luminous and support 
winds with higher velocity and mass loss rates (Crowther et al. l2002> . 

* Two other dust WC9d stars - WR69 & 103 - have also been sug- 
gested to be short period binaries on the basis of periodic photometric 
variability. For the former star Marchenko et al. ( 1998 i report a 2.29 day 



J. S. Clark et al.: The WC9d binary W239 



7 



pact dust forming WC binary in the Galaxy - W239 potentially 
presents an ideal testbed for such a prediction. 

Hydrodynamical simulations of the wind collision zones by 
Pittard ( 12009 ) indicate that coherent large scale Archimedean 
spirals may form in massive short period systems, although 
their properties are highly dependent on the binary separa- 
tion. Specifically, the large scale wind structure of the most 
compact binary considered by Pittard (2009, Porh ~3 days) 
quickly fragmented, leading to a highly clumpy, inhomogeneous 
outflow. Could dust formation occur in such an environment? 
Unfortunately, ab initio simulations of such a process are not 
yet available and are certainly premature for W239. However, 
radiative transfer models of the wider {Porh ~241.5 day) binary 
WR104 of Harries et al. (I2004I I suggest that self shadowing from 
the dense inner regions of wind spiral strongly influences the 
temperature structure and hence location of dust condensation; it 
is unclear whether such a process would also operate in the frag- 
mented wind spiral anticipated for W239, which in turn might 
be expected to affect dust production if the density at which the 
wind cooled to below the dust sublimation temperature was too 
low to permit dust grain growth. 

In this regard we note with interest that a number of WC-i-O 
star CWBs with comparable periods are not dust producer^ al- 
though these are all of earlier spectra type (WC4-7). One might 
speculate that the production of dust in W239 is a result of 
the softer radiation field of the WC9 primary when compared 
to these hotter stars, increasing the likelihood of survival of 
dust grains in this system, although the recent discovery of dust 
asociated with the LMC system HD 36402 (WC4(-hO?)-h08 I; 
Williams 1201 11 1 potentially presents a counterarguement to this 
hypothesis. 

3.3. W239 in an evolutionary context 
3.3.1. Possible binary patliways 

With an orbital period of 5.05 days, W239 is one of the most 
cornpact WC binary system to have been unambiguously identi- 
fiecQ although without dynamical mass estimates for both com- 
ponents we may not yet produce a tailored evolutionary his- 
tory for it. However following the arguements of Petrovic et al. 
( I2005al l. Cantiello et al. (|20^ and De Mink et al. ( |2009] | we 
may conclude that it must have been born in a particularly com- 



periodicity from Hipparcos data but are unable to detemine its origin 
(binarity or NRP), although earlier observations by Balona et al. 1) 1989) 
failed to identify this. Conversely, both the latter authors and Moffat et 
al. ( il986t report a ~1.75 day photometric period for WR 103 which 
was not found in the Hipparcos dataset. Given the potentially transient 
nature of the periodicity in both stars and the lack of a confirmatory 
RV curve we regard the adoption of these periods as orbital in origin as 
premature. 

5 WR 9 (WC7-I-05-8, P„,t, -14.3 day ; Niem ela [T995l l, WR 42 
(WC7+07 V, Po,b -7.9 day; Hill et al. |2000] |. WR79 (WC5+06, 
Porb -8.9 day; Hill et al. 2000 ) and the L MC s ystem HD 36521 
(WC4+06 V-III), PoA - 1 .92 day; Moffat et al. [T990l) 

* In addition to those stars given in footnote 5, van der Hucht ( 1200 It 
lists 3 further WC stars with shorter periods - WR 14, 50 & 86 - in 
all cases solely determined from photometry with no corrobarative RV 
data. As the author cautions, the dataset present by Shylaja ( 11990b for 
WR 14 is insufficient to derive a period and so its identification as a 
binary is highly uncertain. Likewise Veen & Wieringa (2000i cast seri- 
ous doubt upon the photometric variability suggested for WR 50 by van 
Genderen et al. t l991 J . Finally, Paardekooper et al. C2002j show that the 
periodicity associated with WR86 is instead attributable to a yS-Cephei 
companion. 



pact configuration, although none of the systems simulated by 
these authors have such a short orbital period when the primary 
becomes a WR star (Poit >6 days), still less once subsequent 
mass loss has led to a WC phase and caused the orbit to widen. 
Nevertheless, the WN7o-i-? Poa -3.51 day eclipsing binary WR 
B (Bonanos 120071) clearly demonstrates that suitable progenitor 
systems for W239 are present within Wd 1 . 

For evolutionary scenarios invoking non-conservative mass 
transfer (e.g. Petrovic et al. l2005al Cantiello et al. 12007) . increas- 
ing the initial orbital period of the progenitor O-i-O binary re- 
sults in a progressively more compact WR+O MS configuration. 
However, above a critical period and after the transfer of ~50% 
of the envelope mass of the primary, the orbit widens and the 
mass transfer rate increases to the point where the secondary also 
expands to fill its Roche lobe (Wellstein et al. '2001). Adopting 
a standard common envelope prescription at this point, unpub- 
Ushed simulations by Wellstein (2000) found that a I6-1-I5M0, 
Porb ~8 day MS binary evolves to a 5-1-20 Mq, Poit ~ 1 .8 day con- 
figuration. Comparable simulations with improved input physics 
are currently being undertaken by Schneider et al. (in prep.), 
and although a higher mass progenitor would be required (e.g. 
3O...4OM0 for a current mass of -10 M©) such a common enve- 
lope pathway appears promising for the production of W239. 

As an alternative to binary evolution driven by mass trans- 
fer we also highlight the Case M(ixing) evolutionary pathway 
of De Mink et al. (2009). In this scenario - applicable to close, 
massive tidally locked binaries - rapid rotation leads to efficient 
rotational mixing and hence chemically homogeneous evolution. 
This results in the stars evolving directly to the WR phase, while 
remaining within their Roche Lobes. As with the above path- 
ways, once the primary becomes a WR star, wind driven mass 
loss causes the initially compact orbit to widen. However, given 
that we anticipate progenitor and cuiTent masses of -40 & 10 M© 
for the WC9 component of W239, we might expect the orbital 
period to have widened from an initial value of 1.5. ..2. 5 days to 
beyond the observed -5.05 days, although confirmation of this 
would require tailored modeling. 

What implications do these scenarios have on the mass of 
the 5.05 day binary companion? Case B evolution leading to 
a common envelope phase leads to merger for low secondary 
masses, suggesting that it must be rather massive. From the un- 
published thesis work of Wellstein ( I2000I I a time averaged mass 
transfer efficiency of -50% was found for the I6-1-I5M0 binary 
considered. If this is applicable to higher mass systems such as 
W239 we estimate that a -40+38 Mq system would result in 
a current -15-1-50 Mq configuration, which we may exclude on 
observational grounds. For lower initial mass ratios one would 
expect the mass transfer efficiency to be reduced (e.g. Petrovic 
et al. 2005a) and so one could anticipate a -40-1-20 M© system 
evolving into a -15-)-25M0 WR-i-MS system, broadly consistent 
with our cuiTent observational data. Weaker constraints on the 
mass of the binary companion are obtained for a case M evolu- 
tion, but are expected to be in the range -1O-38M0, although 
masses in the upper reaches of this range are again likely to be 
excluded on observational grounds. Nevertheless, while tailored 
numerical modeling is required to confirm these estimate, it ap- 
pears highly likely that the binary companion of W239 is mas- 
sive enough to generate a wind collision zone and support dust 
production. 

3.3.2. Evolutionary pathways in Wdl 

Even if we are unable to provide a unique evolutionary scheme 
for W239, such pathways are of particular interest given that 
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Table 5. Summary of dusty WC binary properties (top panel) and remaining X-ray detected binaries (lower panel). 



WR 


Spectral Type 


Porb 


hx Period determinations 


Reference 








(10^3 ergs-') 





77n/W239 WC9vd+? 



5.053+0.002 day 



19 WC4pd + 09.6 10.1 yr 

48a WC8ed + ? ~1 yr + >30 yr? 



65 


WC9d+? 


'few years' 


70 


WC9vd + BO I 


1045±60 day 


98a 


WC8-9va + .' 


564±8 day 


104 


WC9pd + B0.5 V 


241.5±0.5 day 


112 


WC9d + ? 


24.8±1.5 yr 


113 


WC8d+08-9IV 


-29.707 day 


118 


WC9d + ? 


-60 day 


125 


WC7ed + 09 III 


>4 yr 


137 


WC7ed + 09 


13.05±0.15yr 


140 


WC7pd +04-5 


7.94±0.03 yr 


HD 36402 


WC4(+0?)+08 I 


3.03269 



0.4 



100 



0.1 



5-12 



RV 

photometric 

photometric (provisional) 

variable X-ray 
photometric 
photometric -I- imaging 

imaging 

imaging (provisional) 
RV 

imaging (provisional) 
photometric (provisional) 
RV + photometric 

RV + photometric 
RV 



This work 

Clark et al. ( l2008t 

Veen et al. ( fT998) 



Williams et al. ("2003] 
Zhekov et al. (.2010 1 



Oskinova & Hamann ( I2008t 
Williams et al. ( [2008) 



Monnier et al. CI 999) 
Williams et al. (20031 



Tuthill et al. (|2008] | 
Crowther ([1997} 



Marchenko et al. (12002) 
Massey & Niemela ( I198H 
Millour et al. ( [2007) 
Williams et al. ([1994) 



Williams et al ([2001) 
Lefevre et al. ( [2005) 



Williams et al ([T990) 
Pollock et al. f2005) 
Zhekov & Skinner ( [2000) 



Moffat et al. ([1990) 
Williams ([20TT) 
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48 



79 



WC8+07.5III 78.53+0.01 day 1.2-1.5 



WC5+O9.5/B0Iab 19.138±0.003 day 1.0 
(WC6+06V) 



WC7+05V 



8.891 1±0.0001 day 0.05 



RV 



RV 



RV 



Schmutz et al. ([T997) 
Schild et al. i2004) 



Hill et al. HW^ 
Oskinova (200^ 



Hill et al. f2000) 
Oskinova ( [2005) 



Column 5 summarises the origin of the orbital period, being derived from Radial Velocity data, periodicities in the photometric lightcurve, the 
motion of a pinwheel nebula (imaging) and/or a variable X-ray flux. In several cases (WR48a, 65, 112, 118 & 125) these are less robust and so 
are labeled as provisional. RV variations are also present in WR70 and are consistent with both short (-50 d) and long (-1200 d) periods. Note 
the discrepant classifications given for WR48 in Hill et al. (|2002jl and Oskinva (|2005J|. Naze 12009) report an X-ray detection of the WC7+07/9 
binary WR93 but no X-ray flux is given. 



Wd 1 appears to host a number of compact WR binary sys- 
tems; W239 being the third following WR B (WN7oH-?) and 
W13 (WNVL-i-OB). To these we might also add, albeit with less 
certainty, W44 (WR L; WN9h) - which showed rapid, dramatic 
line profile variability for which binarity is the most compelling 
explanation (Clark et al. 12010) - as well as the progenitor of the 
magnetar CXOU J 1647 10.2-4552 16 (Ritchie et al. i20T0a) . 

However, Wdl 1 also contains a significant population 
of cool supergiants - 10 Yellow Hypergiants (YHGs) & Red 
Supergiants (RSGs) - and a further 3 cool B hypergiants tran- 
siting to/from such a phase (Clark et al. 120051 1201 Oi l. Due to 



their intrinsic luminosity, these stars have such a large physical 
extent (e.g. ~2000 R© for the RSG W26 assuming an intrinsic 
luminosity of log(L/Lo)~5.8) that any putative companion must 
be in such a wide orbit that it could not influence their evolution. 



Consequently, given the apparent coevality of Wd 1 , the twin 
hypergiant and shoit period WR binary populations provide 
compelling evidence for the presence of two distinct evolu- 
tionary pathways for ~3O-4OM0 stars. The first, experienced 
by single stars and long period binaries is predicted by current 
evolutionary theory (e.g. Meynet & Maeder i2003j and results in 
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a post-Main Sequence red loop across the H-R diagram: 

Single channel: 06-7 V ^ 08-9 III ^ 09-B1 la ^ B5-9 
Ia+/YHG ^ RSG ^ B5-9 Ia+/YHG/LBV ^ WN ^ WCQ -> 
SN + production of black hole 

The second, experienced by short period binaries avoids the 
red loop due to binary driven mass transfer: 

Binary channel: 06-7 V + OB V ^ WNLh + OB III-V 
(e.g. Wl^l) ^ WNo + OB III-V (e.g. WR B) ^ WC + OB 
III-V (e.g. WR F) — > SN + production of neutron star 

where one would predict that the wider the initial binary sep- 
aration, the further through a red loop the primary travels before 
the onset of mass transfer curtails it. 

The pre-SN evolutionary pathway followed also plays a sig- 
nificant role in the nature of the post-SN relativistic remnant. 
If binary driven mass loss removes the H-rich mantle in - or 
shortly after - the main sequence, the He-rich core is exposed for 
a longer period of time than in single star evolution, allowing the 
enhanced WR mass loss rates to act for longer This results in a 
lower pre-SN core mass and hence the likely formation of a neu- 
tron star rather than a black hole (Fryer et al. 12002) as appears 
to be the case in Wdl (Clark et al. l2U08l Ritchie et al. I2010al l. 
In contrast, a single star will retain its H-rich mantle for longer, 
be subject to WR mass loss rates for a shorter period of time and 
will likely yield a Black Hole in the mass range considered here 
(e.g. Wellstein & Langer [T999] l. 

Finally case M evolution, which occurs in very compact sys- 
tems (1.5. ..2. 5 days; De Mink et al. 2009) and proceeds without 
significant binary driven mass-loss is thought to lead to rather 
massive pre-SNe core masses and the consequent formation of a 
high mass black hole. 

3.3.3. Wider implications 

It therefore appears that in addition to initial mass, metallicity 
and rotational velocity, binarity plays a key role in the evolu- 
tion of massive stars. If correct, what are the consequences of 
such an hypothesis? Firstly it may help to resolve the difficulties 
in reconciling the current post-MS stellar population of Wd 1 
to evolutionary predictions. Crowther et al. (2006ai highlight 
this problem - specifically that the ratio of WC/WN and WN(H- 
rich)/WN(H-poor) WRs are discrepant - and indeed advance bi- 
narity as a possibly ammeliorative factor 

Secondly, the operation of two parallel evolutionary chan- 
nels would suggest that extrapolating the cluster binary fraction 
from the Wolf-Rayet population may be in error If this were the 
case then one might expect the binary fraction of the unevolved 
O Main Sequence/Giants to be lower than that inferred for the 
WRs; our full multiobject spectroscopic dataset (Ritchie et al. 
12009 ) will allow us to test this assertion. 

Beyond Wd 1, if proven, this hypothesis will have signifi- 
cant impact on the determination of ages and integrated masses 

' We note that the WC8 star WR K appears to be single, being 
amongst the faintest of the WRs within Wdl and demonstrating no ev- 
idence of RV binary reflex motion, an IR excess from dust emission 
or X-ray emission. Hence it appears that single stars with progenitor 
masses ~40 Mq (Ritchie et al. 2010a I can also evolve through a WC 
phase prior to SNe. 

* With a period of ~9.27 days W13 will not evolve into a comparable 
configuration as WR B & F since stellar wind losses will drive it to 
longer periods. 



of star forming regions. As an example we consider the star for- 
mation region recently detected at the base of the Scutum-Crux 
arm and delineated by at least four red supergiant clusters (Figer 
et al. 2006; Davies et al. l2007l Clark et al. 2009, Negueruela et 
al.|2010a). Given that extinction has prevented the identification 
of the Main-Sequence of each cluster, their masses have been de- 
termined via comparison of the number of RSGs to evolutionary 
predictions. If a proportion of massive stars instead avoid this 
phase via close binary evolution then integrated cluster masses 
and consequently star formation rates may be (substantially) un- 
derestimated (Davies et al. 2009). However, as a corollary one 
might also expect a population of (underluminous) WR stars 
produced via the binary channel, which have yet to be identified. 

Indeed, the significance of a close binary pathway in mas- 
sive stellar evolution has already been highlighted by a number 
of authors. Brinchmann et al (12008) found that binary evoluton 
is necessary to replicate the observational properties of Wolf- 
Rayet galaxies at low metalicities. Moreover Smith et al. (1201 Oi l 
argue that the observed fraction of Type Ibc and lib SNe re- 
quires a large proportion of their progenitors to have evolved via 
close binary evolution and mass stripping via Roche Lobe over- 
flow, while Cantiello et al. (2007 ) demonstrate the role such a 
channel will play in the production of GRBs in lower metalicity 
environments. Finally, as implied by the discussion above, the 
prevalence of this pathway will also profoundly influence the 
production rate and properties of both intermediate- and high- 
mass X-ray binaries (e.g. Kobulnicky & Fryer 2007 1, favouring 
the production of neutron star rather than black hole accretors. 

4. Conclusions 

We have analysed multi-epoch spectroscopy and photometry 
of the WC9d star W239, which reveal it to be a binary with 
a 5.05 day period. Non-LTE modeling indicates that the pa- 
rameters of the primary are fully comparable to those of other 
Galactic WC9 stars. Inspection of the spectral energy distri- 
bution of W239 indicates the presence of hot (-1300 K) dust 
which, from the long term near-IR lightcurve appears to be 
present at all epochs. These observations reveal substantial vari- 
ability, with a significant near-IR 'flare' between 2004-6 that is 
best understood as a transient increase in the dust production 
rate resulting from the periastron passage of a massive tertiary 
companion in a wide, eccentric orbit. 

The combination of the current stellar census of Wdl and 
the requirement that they support strong stellar winds constrains 
both secondary and tertiary components to be stars of spectral 
type >07-8 V or 08-9 III. Such classifications are consistent 
with the assumptions made during our non-LTE model atmo- 
sphere analysis, while out tentative identification of a Pal2 pho- 
tospheric line suggests the presence of at least one 07-9 V-III 
star within the W239 system. 

W239 has one of the shortest orbital periods yet determined 
for any WC star This finding leads to several important conclu- 
sions. Firstly, it suggests that dust production in WC stars can 
occur over all binary separations and is not restricted to long pe- 
riod systems, fully consistent with the hypothesis that the persis- 
tent dust producing systems are all short period binaries. Future 
RV observations of the complete dusty WC population within 
Wd 1 will test this scenario. 

Secondly, W239 is the third short period WR binary identi- 
fied within Wd 1 . The compact nature of these systems implies 
that they must have evolved from massive, short period progen- 
itor systems likely via non-conservative mass transfer while still 
on the main sequence, or shortly thereafter In doing so they 
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would have avoided the post-main sequence red-loop that results 
in the formation of the cool YHGs and RSGs that are present 
within Wd 1. We therefore conclude that there is compelUng ev- 
idence for the presence of two parallel evolutionary pathways for 
massive ~3O-4OM0 stars within Wd 1, depending on their binary 
properties (specifically their orbital separation). In combination 
with our current RV survey on OB stars, an extention to include 
all remaining WRs with Wd 1 scheduled for summer 2011 will 
enable us to determine the relative weighting of each channel. 
Consequently, the co-eval stellar population within Wd 1 may 
prove critical in determining whether binarity plays as signifi- 
cant a role in the evolution of massive stars as their initial mass, 
chemical composition and rotational velocity are thought to. 

If this is the case, it will have important consequences for 
the nature of both GRB and SNe progenitors as well as the sub- 
sequent production of relativistic objects - such as the Wd 1 
magnetar - and X-ray binaries. Moreover, this evolutionary bi- 
furcation will need to be incorporated into population synthesis 
codes used to infer the properties of stellar aggregates both in our 
own and external galaxies; for example, mass determinations for 
the RSG-rich clusters at the base of the Scutum-Crux arm may 
currently be significantly underestimated due to this omission. 
Finally, one might speculate that if binaries are predominantly 
formed via dynamical capture in stellar aggregates, rather than 
being primordial, then massive stellar evolution may be partially 
dependent on envirormient, with denser regions leading to higher 
capture probabilities and hence enhanced evolution via the close 
binary channel. 

Acknowledgements. JSC acknowledges support from an RCUK fellowship. 
This research is partially supported by the Spanish Ministerio de Ciencia e 
Innovacion (MICINN) under grants AYA20()8-06166-C03^03,AYA2010-21697- 
C05-()5 and CSD2006-70 and based upon observations made at the Pico dos Dias 
Observatory (LNA/MCT). We wish to thank the referee, Peredur Williams, for 
his comments which greatly improved the paper and Julian Pittard for his valu- 
able input into understanding the likely properties of the wind collision zone in 
W239. 



References 

Allen, D. A., Swings, J. P, Harvey, R M., 1972, A&A, 20, 333 

Balona, L. A., Egan, J., Marang, R, MNRAS, 1989, 103 

Bonanos, A. Z., 2007, AJ, 133, 2696 

Bonnell, I. A., Bate, M. R., 2005, MNRAS, 362, 915 

Bosch, G., Terlevich, E., Terlevich, R., 2009, AJ, 137, 3437 

Brinchmann, J., Kunth, D., Durret, E, 2008, A&A, 485, 657 

Cantiello, M., Yoon, S.-C, Langer, N., Livio, M., 2007, A&A, 465, L29 

Chene, A.-N., St-Louis, N., 2010, ApJ, 716, 929 

Clark, J. S., Negueruela, I., 2002, A&A, 396, L25 

Clark, J. S., Goodwin, S. P, Crowther, P A., et al., 2002, A&A, 392, 909 
Clark, J. S., Negueruela, I., Crowther, P A., Goodwin, S. R, 2005, A&A, 434, 
949 

Clark, J. S., Muno, M. R, Negueruela, I., et al., 2008, A&A, 477, 147 
Clark, J. S., Negueruela, I., Davies, B., et al., 2009, A&A, 498, 109 
Clark, J. S., Ritchie, B. W., Negueruela, I., 2010, A&A, 514, 87 
Cox, A. N., 2001, Astrophysical Quantities 
Crowther, P A., 1997, MNRAS, 290, L59 

Crowther, P A., Dessaart, L., HilUer, D. J., Abbott, J. B., FuUerton, A. W., 2002, 

A&A, 392, 653 

Crowther, P A., Hadfield, L. J., Clark J. S., Negueruela, I., Vacca, W. D., 2006a, 

MNRAS, 372, 1407 
Crowther, R A., Morris, P W., Smith, J. D., 2006b, ApJ, 636, 1033 
Crowther, R A., Barnard, R., Carpano, S., et al., 2010, MNRAS, 403, L41 
Davies, B., Figer, D. E, Kudritzki, R.^R, et al., 2007, ApJ, 671, 781 
Davies, B., Origlia, L., Kudritzki, R.-R, et al, 2009, ApJ, 696, 2014 
De Marco, O., Crowther, P A., Schmutz, W., et al., 2000, A&A, 358, 187 
De Mink, S. E., Cantiello, M., Langer, N., et al., 2009, A&A, 497, 243 
Dougherty, S.M., Clark, J. S., Negueruela, I., Johnson, T., Chapman, J. M., 2010, 

A&A, 511, A58 

Eiger, D. E, MacKenty, J. W., Robberto, M., et al., 2006, ApJ, 643, 1166 
Eryer, C. L., Heger, A., Langer, N., Wellstein, S., 2002, ApJ, 578, 335 



Harries, T. J., Monnier, J. D., Symington, N. H., Kurosawa, R., 2004, MNRAS, 

350, 565 

Hill, G. M., MoflFat, A. F. J., St-Louis, N., Bartzakos, R, 2000, MNRAS, 318, 

402 

Hill, G. M., Moffat, A. R J., St-Louis, N., 2002, MNRAS, 335, 1069 
Hill, G. M., Moffat, A. E J., St-Louis, N., Bartzakos, R, 2000, MNRAS, 318, 
402 

Hillier, D. J., Miller, D. L., 1998, ApJ, 496, 407 
Hillier, D. J., Miller, D. L., 1999, ApJ, 519, 354 
van der Hucht, K. A., 2001, NewAR, 45, 135 

van der Hucht, K. A., Williams, R M., Morris, P W., 2001, ESASP, 460, 273 
Kaper, L., van der Meer, A., Najarro, E, 2006, A&A, 457, 595 
Kiminki, D. C, Kobulnicky, H. A., Kinemuchi, K., et al., 2007, ApJ, 664, 1102 
Kobulnicky, H. A., Eryer, C. L., 2007, ApJ, 670, 747 
Lefevre, L., Marchenko, S. V., Lepine, S., et al., 2005, MNRAS, 360, 141 
Marchenko, S. V., Moffat, A. E J., van der Hucht, K. A., et al., 1998, A&A, 331, 
1022 

Marchenko S. V., Moffat, A. F. J., Vacca, W. D., Cote, S., Doyon, R., 2002, ApJ, 
565, 59 

Marchenko, S. V., Moffat, A. F J., 2007, ASRC, 367, 213 

Martins, F, Schaerer, D., Hillier, D., 2005, A&A, 436, 1049 

Martins, F, Plez, B., 2006, A&A, 457, 637 

Massey, P, Niemela, V. S., 1981, ApJ, 245, 195 

Mengel, S., Tacconi-Garman, L. E., 2009, Ap&SS, 324, 321 

Meynet, G., Maeder, A., 2003, A&A, 404, 975 

Millour, E, Petrov, R. G., Chesneau, O., et al., 2007, A&A, 464, 107 

de Mink, S. E., Cantiello, M., Langer, N., et al., 2009, A&A, 497, 243 

Moffat, A. F. J., Lamontagne, R., Cerruti, M., 1986, PASR 98, 1170 

Moffat, A. F J., Niemela, V. S., Marraco, H. G., 1990, ApJ, 348, 232 

Monnier, J., Tuthill, R, Danchi, W., 1999, ApJ, 525, L97 

Morel, T., St-Louis, N., Marchenko, V., 1997, 482, 470 

Naze, ¥., 2009, A&A, 506, 1055 

Negueruela, I., Gonzalez-Fernandez, C, Marco, A., Clark, J. S., Martfnex- 

Nufiez, S., 2010a, A&A, 511, 840 
Negueruela, I., Clark, J. S., Ritchie, B. W., 2010b, A&A, 516, A78 
Niemela, V. S., 1995, lAUS, 163, 223 

Oskinova, L. M., Ignace, R., Hamann, W.-R., Pollock, A. M. T., Brown, J. C, 

2003, A&A, 402, 755 
0.skinova, L. M., 2005, MNRAS, 361, 679 
Oskinova, L. M., Hamann, W.-R., 2008, MNRAS, 390, L78 
Owocki, S. P, Gayley, K. G., 1995, ApJ, 454, 1250 

Paardekooper, S. J., Veen, P. M., van Genderen, A. M., van der Hucht, K. A., 

2002, A&A, 384, 1012 
Raczynski, B., 1967, Acta Astron., 17, 355 

Rasquani, L., Avila, G., Blecha, A., et al., 2002, The Messenger, 110, 1 
Retrovic, J., Langer, N., van der Hucht, K. A., 2005a, A&A, 435, 1013 
Retrovic, J., Langer, N., Yoon, S.-C, Heger, A., 2005b, A&A, 435, 247 
Riatti, A. E., Bica, E., Claria, J. J., 1998, A&AS, 127, 423 
Rittard, J. M., 2009, MNRAS, 396, 1743 
Rodsiadlowski, P, Joss, R C, Hsu, J. J. L., 1992, ApJ, 391, 246 
Pollock, A. M. T., Corcoran, M. F., Stevens, I. R., Williams, P. M., 2005, ApJ, 
629, 482 

Prilutskii, O. F, Usov, V. V., 1976, Sov. Astr, 20, 2 

Ritchie, B. W., Clark, J. S., Negueruela, I., Crowther, P A., 2009, A&A, 507, 
1585 

Ritchie, B. W., Clark, J. S., Negueruela, I., Langer, N., 2010a, A&A, 520, A48 

Ritchie, B. W., Clark, J. S., Negueruela, I., 2010b, BSRSL 80, 628 

Sana, H., Gosset, E., Naze, Y, Rauw, G., Linder, N., 2008, MNRAS, 386, 447 

Schaerer, D., Maeder, A., 1992, A&A, 263, 129 

Schild, H., Gudel, M., Mewe, R., et al., 2004, A&A, 422, 177 

Schmutz, W., Schweickhardt, J., Stahl, O., et al., 1997, A&A, 328, 219 

Shylaja, B. S., 1990, Ap&SS, 164, 63 

Silverman, J. M., Eilippenko, A. V, 2008, ApJ, 678, L17 

Skinner, S. L., Simmons, A. E., Zhekov, S. A., et al., 2006, ApJ, 639, L35 

Smith, N., Li, W., FiUppenko, A. V., Chomock, R., 2010, MNRAS, in press. 

Stevens, I. R., Howarth, I. D., MNRAS, 302, 549 

Tuthill, P G., Monnier, J. D., Danchi, W. C, 1999, Nature, 398, 487 

Tuthill, P G., Monnier, J. D., Lawrance, N., et al., 2008, ApJ, 675, 698 

Usov, V. v., 1991, MNRAS, 252, 49 

van Genderen, A. M., Verheijen, M. A. W., van der Hucht, K. A., et al., 1991, 
lAUS, 143, 129 

Veen, R M., van der Hucht, K. A., Williams, R M., et al., 1998, A&A, 339, L45 
Veen, P M., Wieringa, M. H., 2000, A&A, 363, 1026 
Wellstein, S., 2000, PhD thesis, Potsdam University 
Wellstein, S., Langer, N., 1999, A&A, 350, 148 
Wellstein, S., Langer, N., Braun, H., 2001, A&A, 369, 939 
Wiffiams, R M., van der Hucht, K. A., Pollock, A. M. T., et al., 1990, MNRAS, 
243, 662 



J. S. Clark et al.: The WC9d binary W239 



11 



WilUams, P. M., van der Hucht, K. A., 1992, in A. S. P. Conf. Series, Vol. 22, 
269 

Williams, P. M., van der Hucht, K. A., Kidger, M. R., Geballe, T. R., Bouchet, 

P, 1994, MNRAS, 266, 247 
Williams, R M., van der Hucht, K. A., 2000, MNRAS, 314, 23 
Williams, R M., Kidger, M. R., van der Hucht, K. A., et al., 2(X)1, MNRAS, 324, 

156 

Williams, R M., van der Hucht, K. A., Morris, P W., Marang, F., 2003, lAUS, 
212, 115 

Williams, P M., 2008, RevMexAA, 33, 71 

WilUams, P. M., Marchenko, S. V., Marston, A. P, et al., 2009, MNRAS, 395, 
1749 

Williams, R M., 201 1, BSRSL, 80, 195 

Zhekov, S. A., Skinner, S. L., 2000, Apl, 538, 808 

Zhekov, S A., Gagne, M., Skinner, S. L., ApJL, 2011, 727, L17 

Zinnecker, H., Yorke, H. W., 2007, ARA&A, 45, 481 



